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We study by analytical methods and large scale simulations a dynamical model for the spreading 
of epidemics in complex networks. In networks with exponentially bounded connectivity we recover 
the usual epidemic behavior with a threshold defining a critical point below which the infection 
prevalence is null. On the contrary, on a wide range of scale-free networks we observe the absence 
of an epidemic threshold and its associated critical behavior. This implies that scale-free networks 
are prone to the spreading and the persistence of infections whatever spreading rate the epidemic 
agents might possess. These results can help understanding computer virus epidemics and other 
spreading phenomena on communication and social networks. 
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I. INTRODUCTION 

Many social, biological, and communication systems 
can be properly described by complex networks whose 
nodes represent individuals or organizations and links 
mimic the interactions among them . Recently, many 
authors have recognized the importance of local cluster- 
ing in complex networks. This implies that some special 
nodes of the network posses a larger probability to de- 
velop connections pointing to other nodes. Particularly 
interesting examples of this kind of behavior are found 
in metabolic networks food webs Q, and, most im- 
portantly, in the Internet and the world- wide- web, where 
the networking properties have been extensively studied 
because of their technological and economical relevance 

Complex networks can be classified in two main 
groups, depending on their connectivity properties. The 
first and most studied one is represented by the expo- 
nential networks, in which the nodes' connectivity dis- 
tribution (the probability P{k) that a node is connected 
to other k nodes) is exponentially bounded [p|-pO|. The 
typical example of an exponential network is the random 
graph model of Erdos and Renyi Q . A network belong- 
ing to this class that has recently attracted a great deal of 
attention is the Watts and Strogatz model (WS) 0-|2), 
which has become the prototypical example of a small- 
world network ||l3|] . A second and very different class of 
graph is identified by the scale-free (SF) networks which 
exhibit a power-law connectivity distribution |p^. 



P(fc) - fc~2-7^ 



(1) 



where the parameter 7 must be larger than zero to ensure 
a finite average connectivity (fe). This kind of distribu- 
tion implies that each node has a statistically significant 
probability of having a very large number of connections 
compared to the average connectivity of the network (fc). 
In particular, we will focus here on the Barabasi and 



Albert model (BA) |14| , which results in a connectivity 
distribution P{k) ^ . 

In view of the wide occurrence of complex networks in 
nature, it becomes a very interesting issue to inspect the 
effect of their complex features on the dynamics of epi- 
demic and disease spreading |15|, and more in general on 
the nonequilibrium phase transitions that usually charac- 
terize these type of phenomena jl^]. It is easy to foresee 
that the characterization and understanding of epidemic 
dynamics on these networks can find immediate applica- 
tions to a wide r ang e of problems, r ang ing from computer 
virus infections |17| ] , epidemiology [ |l8[ , and the spreading 
of polluting agents jl^ . 

In this paper, we shall study the susceptible-infected- 
susceptible (SIS) model jl^ on complex networks. We 
study analytically the prevalence and persistence of in- 
fected individuals in exponential and SF networks by us- 
ing a single-site approximation that takes into account 
the inhomogeneity due to the connectivity distribution. 
We find that exponential networks show, as expected, an 
epidemic threshold (critical point) separating an infected 
from a healthy phase. The density of infected nodes de- 
creases to zero at the threshold with the linear behavior 
typical of a mean-field (MF) critical point jl6|. The SF 
networks, on the other hand, show a very different and 
surprising behavior. For < 7 < 1 the model does not 
show an epidemic threshold and the infection can always 
pervade the whole system. In the region 1 < 7 < 2, the 
model shows an epidemic threshold that is approached, 
however, with a vanishing slope; i.e. in the absence of 
critical fluctuations. Only for 7 > 2 we recover again the 
usual critical behavior at the threshold. In these systems, 
because of the nonlocal connectivity, single site approx- 
imation predictions are expected to correctly depict the 
model's behavior. In order to test our predictions, we 
perform large scale numerical simulations on both expo- 
nential and SF networks. Numerical results are in perfect 
agreement with the analytical predictions and confirm 
the overall picture for the SIS model on complex networks 
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given by the theoretical analysis. The striking absence of 
an epidemic threshold on SF networks, a characteristic 
element in mathematical epidemiology, radically changes 
many of the conclusions drawn in classic epidemic model- 
ing. The present results could be relevant also in the field 
of absorbing-phase transitions and catalytic reactions in 
which the spatial interaction of the reactants can be mod- 
eled by a complex network [l^ . 

The paper is organized as follows. In Sec. II we intro- 
duce the SIS model in a general context. Sec. Ill is de- 
voted to the analysis of exponentially bounded networks, 
exemplified by the WS model. In Sec. IV we analyze 
the scale- free BA model, with connectivity P{k) ~ k~^. 
Sec. V extends the analytical approach applied to the 
BA model to generalized SF networks, with connectivity 
distribution P{k) - k-^"'' , 7 > 0. Finally, in Sec. VI 
we draw our conclusions and perspectives. 

II. THE SIS MODEL 

To address the effect of the topology of complex net- 
works in epidemic spreading we shall study the stan- 
dard susceptible-infected-susceptible (SIS) epidemiolog- 
ical model . Each node of the network represents an 
individual and each link is a connection along which the 
infection can spread to other individuals. The SIS model 
relies on a coarse grained description of the individuals in 
the population. Within this description, individuals can 
only exist in two discrete states, namely, susceptible, or 
"healthy" , and infected. These states completely neglect 
the details of the infection mechanism within each indi- 
vidual. The disease transmission is also described in an 
effective way. At each time step, each susceptible node 
is infected with probability v if it is connected to one or 
more infected nodes. At the same time, infected nodes 
are cured and become again susceptible with probability 
i5, defining an effective spreading rate A = 1^/6. (Without 
lack of generality, we set (5 = 1.) Individuals run stochas- 
tically through the cycle susceptible infected sus- 
ceptible, hence the name of the model. The updating can 
be performed with both parallel or sequential dynamics 
Ip^ . The SIS model does not take into account the pos- 
sibility of individuals removal due to death or acquired 
immunization ||l8|| . It is mainly used as a paradigmatic 
model for the study of infectious disease that leads to an 
endemic state with a stationary and constant value for 
the prevalence of infected individuals, i.e., the degree to 
which the infection is widespread in the population. 

The topology of the network that specifies the inter- 
actions among individuals is of primary importance in 
determining many of the model's features. In standard 
topologies the most significant result is the general pre- 
diction of a nonzero epidemic threshold Ac [0. If the 
value of A is above the threshold, A > Ac, the infection 
spreads and becomes persistent in time. Below it, A < Ac, 
the infection dies out exponentially fast. In both sides of 



the phase diagram it is possible to study the behavior in 
time of interesting dynamical magnitudes of epidemics, 
such as the time survival probability and the relaxation 
to the healthy state or the stationary endemic state. In 
the latter case, if we start from a localized seed we can 
study the epidemic outbreak preceding the endemic sta- 
bilization. From this general picture, it is natural to con- 
sider the epidemic threshold as completely equivalent to 
a critical point in a nonequilibrium phase transition . 
In this case, the critical point separates an active phase 
with a stationary density of infected nodes (an endemic 
state) from an absorbing phase with only healthy nodes 
and null activity. In particular, it is easy to recognize 
that the SIS model is a generalization of the contact pro- 
cess (CP) model, that has been extensively studied in the 
context of absorbing-state phase transitions [ p^ . 

In order to obtain an analytical understanding of the 
SIS model behavior on complex networks, we can apply 
a single site dynamical MF approach, that we expect to 
recover exactly the model's behavior due to the nonlocal 
connectivity of these graphs. Let us consider separately 
the case of the exponentially bounded and SF networks. 

III. EXPONENTIAL NETWORKS: THE 
WATTS-STROGATZ MODEL 

The class of exponential networks refers to random 
graph models which produce a connectivity distribution 
P{k) peaked at an average value (fc) and decaying expo- 
nentially fast for k ^ (k) and fc ^ (fc). Typical examples 
of such a network are the random graph model [|9| an d the 
small- world model of Watts and Strogatz (WS) ||l^. The 
latter has recently been the object of several studies as 
a good candidate for the modeling of many realistic sit- 
uations in the context of social and natural networks. In 
particular, the WS model shows the "small-world" prop- 
erty common in random graphs [T^ ]; i.e., the diameter 
of the graph — the shortest chain of links connecting any 
two vertices — increases very slowly, in general logarithmi- 
cally with the network size On the other hand, the 
WS model has also a local structure (clustering property) 
that is not found in random graphs with finite connectiv- 
ity The WS graph is defined as follows 0,|l2|: 
The starting point is a ring with N nodes, in which 
each node is symmetrically connected with its 2K near- 
est neighbors. Then, for every node each link connected 
to a clockwise neighbor is rewired to a randomly chosen 
node with probability p, and preserved with probability 
1 — p. This procedure generates a random graph with a 
connectivity distributed exponentially for large k [|l^,|l2| , 
and an average connectivity (fc) — 2K. The graphs has 
small-world properties and a non-trivial "clustering co- 
efficient"; i.e., neighboring nodes share many common 
neighbors ||l^,|l^. The richness of this model has stim- 
ulated an intense activity aimed at understanding the 
network's properties upon changing p and the network 
size N |l0|-p|,p0|j2l| . At the same time, the behavior 
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of physical models on WS graphs has been investi gated, 
including epidemiological percolation models |ll5| , po| , ^ 
and models with epidemic cycles p3| . 

Here we focus on the WS model with p — 1; it is worth 
noticing that even in this extreme case the network re- 
tains some memory of the generating procedure. The 
network, in fact, is not locally equivalent to a random 
graph in that each node has at least K neighbors. Since 
the fluctuations in the connectivity are very small in the 
WS graph, due to its exponential distribution, we can 
approach the analytical study of the SIS model by con- 
sidering a single MF reaction equation for the density of 
infected nodes p{t): 



dtp{t) = -p{t) + A (fc) p{t) [1 - p{t)] + h.o. 



(2) 



The MF character of this equation stems from the fact 
that we have neglected the density correlations among 
the different nodes, independently of their respective con- 
nectivities. In Eq. (H) we have ignored all higher order 
corrections in p{t), since we are interested in the on- 
set of the infection close to the phase transition; i.e., at 
p{t) <^ 1. The first term on the r.h.s. in Eq. (H) consid- 
ers infected nodes become healthy with unit rate. The 
second term represents the average density of newly in- 
fected nodes generated by each active node. This is pro- 
portional to the infection spreading rate A, the number 
of links emanating from each node, and the probability 
that a given link points to a healthy node, [1 — p{t)]. In 
these models, connectivity has only exponentially small 
fluctuations ((fc^) ~ (fc)) and as a first approximation 
we have considered that each node has the same number 
of links, fc ~ (k). This is equivalent to an homogeneity 
assumption for the system's connectivity. After impos- 
ing the stationarity condition dtp(t) = 0, we obtain the 
equation 



p[-l + A(fc)(l-p)] =0 
for the steady state density p of infected nodes. 



equation defines an epidemic threshold Ac 
yields: 



(fc)- 



p = 

P X- Xc 



if A < Ac 
if A > Ac 



(3) 

This 
, and 



(4a) 
(4b) 



In analogy with critical phenomena, we can consider p 
as the order parameter of a phase transition and A as 
the tuning parameter, recovering a MF critical behavior 
p^ . It is possible to refine the above calculations by 
introducing connectivity fluctuations (as it will be done 
later for SF networks, see Sec. IV). However, the results 
are qualitatively and quantitatively the same as far as 
we are only concerned with the model's behavior close to 
the threshold. 

In order to compare with the analytical prediction we 
have performed large scale simulations of the SIS model 
in the WS network with p = 1. Simulations were im- 
plemented on graphs with number of nodes ranging from 
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FIG. 1. Density of infected nodes p as a function of A in 
the WS network (full line) and the BA network (dashed line). 
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FIG. 2. Log-log plot of density of infected node p as a func- 
tion of A - Ac in WS network, with Ac = 0.1643 ± 0.01. The 
full line is a fit to the form p ~ (A — Ac)^, with an exponent 
fi = 0.97 ±0.04. 



N = 10^ to = 3 X 10^, analyzing the stationary prop- 
erties of the density of infected nodes p; i.e. the infection 
prevalence. Initially we infect half of the nodes in the net- 
work, and iterate the rules of the SIS model with parallel 
updating. In the active phase, after an initial transient 
regime, the systems stabilize in a steady state with a con- 
stant average density of infected nodes. The prevalence 
is computed averaging over at least 100 different starting 
conflgurations, performed on at least 10 different real- 
ization of the random networks. In our simulations we 
consider the WS network with parameter AT = 3, which 
corresponds to an average connectivity (k) = 6. 

As shown in Figs. |l| and ^, the SIS model on a WS 
graph exhibits an epidemic threshold Ac = 0.1643 ± 0.01 
that is approached with linear behavior by p. The value 
of the threshold is in good agreement with the MF pre- 
dictions Ac — 1/2K = 0.1666, as well as the density of 
infected nodes behavior. In Fig. ^ we plot p as a func- 
tion of A — Ac in log-log scale. A linear fit to the form 
p ^ {X — Xc)^ provides an exponent f3 — 0.97 ± 0.04, 
in good agreement with the analytical finding of the 
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FIG. 3. Density of infected nodes p{t) as a function of time 
in supercritical spreading experiments in the WS network. 
Network size N — 1.5 x 10®. Spreading rates range from 
A - Ac = 0.002 to 0.0007 (top to bottom). 
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FIG. 4. Density of infected nodes p{t) as a function of 
time in subcritical spreading experiments in the WS net- 
work. Network size TV = 3 x 10®. Spreading rates range 
from Ac - A = 0.005 to 0.03 (right to left). 



Eq. (^. 

To complete our study of the SIS model in the WS 
network, we have also analyzed the epidemic spreading 
properties, computed by considering the time evolution 
of infections starting from a very small concentration of 
infected nodes. In Fig. || we plot the evolution of the in- 
fected nodes density as a function of time for epidemics 
in the supercritical regime (A > Ac) which start from 
a single infected node. Each curve represents the aver- 
age over several spreading events with the same A. We 
clearly notice a spreading growth faster than any power- 
law, in agreement with Eq. (|^) which predicts an ex- 
ponential saturation to the endemic steady state. In 
the subcritical regime (A < Ac), by introducing a small 
perturbation to the stationary state p — 0, and keep- 
ing only first order terms in Eq. (^, we obtain that 
the infection decays following the exponential relaxation 
dtp{t) = — (k) (Ac — X)p{t). This equation introduces a 
characteristic relaxation time 
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FIG. 5. Surviving probability Ps{t) as a function of time 
in subcritical spreading experiments in the WS network. 
Network size A'^ = 3 x 10®. Spreading rates range from 
Ac — A = 0.005 to 0.03 (right to left). Inset: Surviving prob- 
ability for a fixed spreading rate Ac — A = 0.005 and network 
sizes N = 3x 10^ 10®, and 3 x 10®. 



(k) (Ac -A), 



(5) 



that diverges at the epidemic threshold. Below the 
threshold, the epidemic outbreak dies within a finite 
time; i.e. it does not reach a stationary endemic state. 
In Fig. ^ we plot the average of p{t) for epidemics start- 
ing with an initial concentration po = 0.01 of infected 
nodes; the Figure shows a clear exponential approach to 
the healthy (absorbing) state as predicted by Eq. (||). 
In the subcritical regime, we can compute also the sur- 
viving probability, Ps{t), defined as the probability that 
an epidemic outbreak survives up to the time t [|l6| . In 
Fig. ^ we plot the survival probability computed from 
simulations starting with a single infected node in a WS 
graph of size TV = 3 x 10^. The survival probability decay 
is obviously governed by the same exponential behavior 
and characteristic time of the density of infected nodes 
as confirmed by numerical simulations. Indeed, below 
the epidemic threshold, the relaxation to the absorbing 
state does not depend on the network size N (see inset in 
Fig. |), and the average lifetime corresponding to each 
spreading rate A can be measured by the slope of the 
exponential tail of Ps{t) and p{t). By plotting as a 
function of Ac — A (see Fig. ph, we recover the analytic 
predictions ; i.e. the linear behavior and the unique char- 
acteristic time for both the density and survival proba- 
bility decay. The slope of the graph, measured by means 
of a least squares fitting, provides a value of 6.3, whereas 
the intercept yields 1.0, in good agreement with the the- 
oretical predictions of Eq. (||), (fc) = 6 and (fc) Ac = 1, 
respectively. 

In summary, numerical and analytical results confirms 
that for WS graphs, the standard epidemiological picture 
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FIG. 6. Inverse relaxation time for the SIS model in the 
WS graph as a function of the spreading rate A, estimated 
from the slope of the exponential decay of the infected nodes 
density p{t) (o) and the survival probabihty Ps{t) (O). 



(often called the deterministic approximation) is qualita- 
tively and quantitatively correct. This result, that is well- 
known for random graphs, holds also in the WS model 
despite the different local structure. 



IV. SCALE-FREE NETWORKS: THE 
BARABASI-ALBERT MODEL 

The Barabasi- Albert (BA) graph was introduced as a 
model of growing network (such as the Internet or the 
world-wide-web) in which the successively added nodes 
establish links with higher probability pointing to already 
highly connected nodes This is a rather intuitive 

phenomenon on the Internet and other social networks, 
in which new individuals tend to develop more easily con- 
nections with individuals which are already well-known 
and widely connected. The BA graph is constructed us- 
ing the following algorithm |l^: We start from a small 
number toq of disconnected nodes; every time step a new 
vertex is added, with m links that are connected to an 
old node i with probability 



n(fc.) 



(6) 



where ki is the connectivity of the z-th node. After iter- 
ating this scheme a sufficient number of times, we obtain 
a network composed by A'' nodes with connectivity dis- 
tribution P{k) ~ and average connectivity (fc) = 2rn 
(in this work we will consider the parameters toq = 5 a-nd 
m — 3). Despite the well-defined average connectivity, 
the scale invariant properties of the network turns out to 
play a major role on the properties of models such as per- 
colation p^ , p5[ , used to mimic the resilience to attacks 
of a network. For this class of graphs, in fact, the ab- 
sence of a characteristic scale for the connectivity makes 
highly connected nodes statistically significant, and in- 
duces strong fluctuations in the connectivity distribution 
which cannot be neglected. In order to take into account 



these fluctuations, we have to relax the homogeneity as- 
sumption used for exponential networks, and consider the 
relative density pk {t) of infected nodes with given connec- 
tivity fc; i.e., the probability that a node with k links is 
infected. The dynamical MF reaction rate equations can 
thus be written as 



dtpkit) = -pk{t) + Afc [1 - pk{t)] Q{p{t)), 



(7) 



where also in this case we have considered a unitary re- 
covery rate and neglected higher order terms {p{t) 1). 
The creation term considers the probability that a node 
with k links is healthy [1 — pu (i)] and gets the infection 
via a connected node. The probability of this last event is 
proportional to the infection rate A, the number of con- 
nections k, and the probability Q{p[t)) that any given 
link points to an infected node. We make the assump- 
tion that O is a function of the total density of infected 
nodes p{t) In the steady (endemic) state, p is just a 
function of A. Thus, the probability Q becomes also an 
implicit function of the spreading rate, and by imposing 
stationarity [dtPk{i) = 0], we obtain 



Pk 



kXejX) 
l + kXe{X)' 



(8) 



This set of equations show that the higher the node con- 
nectivity, the higher the probability to be in an infected 
state. This inhomogeneity must be taken into account in 
the computation of 6(A). Indeed, the probability that a 
link points to a node with s links is proportional to sP{s). 
In other words, a randomly chosen link is more likely to 
be connected to an infected node with high connectivity, 
yielding the relation 



e(A) = ^ 



kP{k)pk 



(9) 



Since pk is on its turn a function of 0(A), we obtain a 
self-consistency equation that allows to find 8(A) and an 
explicit form for Eq. (^). Finally, we can evaluate the 
order parameter (persistence) p using the relation 



P 



(10) 



In order to perform an explicit calculation for the BA 
model, we use a continuous k approximation that al- 
lows the practical substitution of series with integrals 
iQ . The full connectivity distribution is given by P{k) = 
2m? /k~^ , where m is the minimum number of connection 
at each node. By noticing that the average connectivity 
is {k) = kP{k)dk = 2m, Eq. (|) gives 



e(A) = mAe(A) 



which yields the solution 



1 



fc3 l + fcAe(A)^ 



(11) 
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FIG. 7. Persistence p as a function of 1/A for BA networks 
of different sizes: = 10^ (+), iV = 5 x 10^ (□), N = 10^ 
(x), iV = 5 X 10'' (o), and iV = 8.5 x 10*^ (O). The linear 
behavior on the semi- logarithmic scale proves the stretched 
exponential behavior predicted for the persistence. The full 
line is a fit to the form p ~ exp(— C/A). 
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(12) 



In order to find the behavior of the density of infected 
nodes we have to solve Eq. (|l^), that reads as 



1 



k 



fc3 l + fcAe(A)' 



(13) 



By substituting the obtained expression for 8(A) and 
solving the integral we find at the lowest order in A 



-l/mA 



(14) 



This result shows the surprising absence of any epi- 
demic threshold or critical point in the model; i.e., Ac — 
0. This can be intuitively understood by noticing that 
for usual lattices and MF models, the higher the node's 
connectivity, the smaller the epidemic threshold. In the 
BA network the unbounded fluctuations in the number 
of links emanating from each node ((fc^) = oo) plays 
the role of an infinite connectivity, annulling thus the 
threshold. This implies that infections can pervade a BA 
network, whatever the infection rate they have. 

The numerical simulations performed on the BA net- 
work confirm the picture extracted from the analytic 
treatment. We consider the SIS model on BA networks 
of size ranging from N = 10'^ to iV = 8.5 x 10^. In Fig. |l| 
we have plotted the epidemic persistence p as a func- 
tion of A in a linear scale. The function p approaches 
smoothly the value A = with vanishing slope. Fig. |^, 
in fact, shows that the infection prevalence in the steady 
state decays with A as p ~ exp(— C/A), where C is a 
constant. The numerical value obtained = 2.5 is 
also in good agreement with the theoretical prediction 
= TO = 3. In order to rule out the presence of fi- 
nite size effect hiding an abrupt transition (the so-called 
smoothing out of critical points |l6[), we have inspected 
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FIG. 8. The density pk, defined as the fraction of nodes 
with connectivity k which are infected, in a BA network of 
size = 5 X 10^ and spreading rates A = 0.1, 0.08 and 0.065 
(bottom to top). The plot recovers the form predicted in 
Eq.d). 
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FIG. 9. Density of infected nodes p{t) as a function of time 
in supercritical spreading experiments in the BA network. 
Network size A'' — 10^. Spreading rates range from A = 0.05 
to 0.065 (bottom to top). 



the behavior of the stationary persistence for network 
sizes varying over three orders of magnitude. The total 
absence of scaling of p and the perfect agreement for any 
size with the analytically predicted exponential behavior 
allows us to definitely confirm the absence of any finite 
epidemic threshold. In Fig. |[ we also provide an illus- 
tration of the behavior of the probability pk that a node 
with given connectivity k is infected. Also in this case 
we found a behavior with k in complete agreement with 
the analytical prediction of Eq. (^) . 

Our numerical study of the spreading dynamical prop- 
erties on the BA network is reported in Figs. ^ and |l^. 
In Fig. ^ we plot the growth of the epidemics starting 
from a single infected node. We observe that the spread- 
ing growth in time has an algebraic form, as opposed to 
the exponential growth typical of mean-field continuous 
phase transitions close to the critical point ||l6| , and the 
behavior of the SIS model in the WS graph (see Fig. |^) . 
The surviving probability Pg (t) for a fixed value of A and 
networks of different size N is reported in Fig HO. In this 
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FIG. 10. Surviving probability Ps{t) as a function of time 
in subcritical spreading experiments in the BA network. 
Spreading rate A = 0.065. Network sizes ranging from 
N = 6.25 X 10^ to TV = 5 X 10^ (bottom to top). 



case, we recover an exponential behavior in time, that 
has its origin in the finite size of the network. In fact, 
for any finite system, the epidemic will eventually die out 
because there is a finite probability that all individuals 
cure the infection at the same time. This probability is 
decreasing with the system size and the lifetime is infinite 
only in the thermodynamic limit N —>■ oo. However, the 
lifetime becomes virtually infinite (the metastable state 
has a lifetime too long for our observation period) for 
enough large sizes that depend upon the spreading rate 
A. This is a well-known feature of the survival probabil- 
ity in finite size absorbing-state systems poised above the 
critical point. In our case, this picture is confirmed by 
numerical simulations that shows that the average life- 
time of the survival probability is increasing with the 
network size for all the values of A. Given the intrinsic 
dynamical nature of scale-free networks, this result could 
possibly have several practical implications in the study 
of epidemic spreading in real growing networks. 

The numerical analysis supports and confirms the an- 
alytical results pointing out the existence of a different 
epidemiological framework for SF networks. The absence 
of an epidemic threshold, a central element in the theory 
of epidemics, opens a different scenario and rationaliza- 
tion for epidemic events in these networks. The danger- 
ous absence of the epidemic threshold, that leaves SF 
networks completely disarmed with respect to epidemic 
outbreaks, is fortunately balanced from a corresponding 
exponentially low prevalence at small spreading rates. In 
addition, the absence of a critical threshold, and the asso- 
ciated diverging response function, makes the increase of 
the endemic prevalence with the spreading rate very slow. 
This new perspective seems to be particularly relevant in 
the rationalization of epidemic data from computer virus 
infections |2^. 



V. GENERALIZED SCALE-FREE NETWORKS 

Recently there has been a burst of activity in the mod- 
eling of SF complex network. The recipe of Barabasi 
and Albert |l^ has been followed by several variations 
and generalizations [p8[- Bl| a nd the revamping of previ- 
ous mathematical works jS^ . All these studies propose 
methods to generate SF networks with variable exponent 
7. The analytical treatment presented in the previous 
section for the SIS model can be easily generalized to SF 
networks with connectivity distribution with 7 > 0. Con- 
sider a generalized SF network with a normalized connec- 
tivity distribution given by 



P{k) = (1 + 7)mi+'^fc-2-7^ 



(15) 



where we are approximating the connectivity k as a, con- 
tinuous variable and assuming m the minimum connec- 
tivity of any node. The average connectivity is thus 



(fc) = / kP{k)dk = ^-^"1. 



(16) 



For any connectivity distribution, the relative density of 
infected nodes is given by Eq. (^. Applying then 
Eq. (^) to compute self-consistently the probability O, 
we obtain 



e(A) = F[1, 7,1 + 7, -(mAe(A))- 



(17) 



where F is the Gauss hypergeometric function p3| . On 
the other hand, the expression for the density p, Eq. (p^), 
yields 



p = F[l,l + 7,2 + 7,-(mAe(A))- 



(18) 



In order to solve Eqs. ( [l7| ) and ( [T^ ) in the limit p ^ 
(which obviously corresponds also to — > 0, we must 
perform a Taylor expansion of the hypergeometric func- 
tion. The expansion for Eq. (n^) has the form p3[ 



i^[l,7,l + 7,-("iAe(A))-i] 

00 

sm(77r) ^ ' I I 



n — 7 



(19) 



where r(a;) is the standard gamma function. An anal- 
ogous expression holds for (|lj). The expansion (|l^) is 
valid for any 7 7^ 1, 2, 3, . . .. Integer values of 7 must be 
analyzed in a case by case basis. (The particular value 
7=1 was considered in the previous Section.) For all 
values of 7, the leading behavior of Eq. (p^) is the same: 



1 + 7 
7 

7 



lAe 



(20) 



The leading behavior in the r.h.s. of Eq. (|T^), on the 
other hand, depends of the particular value of 7: 
(i) < 7 < 1: In this case, one has 



7 



from which we obtain 
77r 



77r 



sin(77r) 



0(A) 



1/(1-7) 



(mA) 



7/(1-7) 



sin(77r) 

Combining this with Eq. (pO[), we obtain; 



(21) 



(22) 



(23) 



Here we have again the total absence of any epidemic 
threshold and the associated critical behavior, as we have 
already shown for the case 7 = 1. In this case, however, 
the relation between p and A is given by a power law with 
exponent /5 = 1/(1 — 7); i.e., /3 > 1. 

(ii) 1 < 7 < 2: In this case, to obtain a nontrivial in- 
formation for 0, we must keep the first two most relevant 
terms in Eq. (nT 



e(A) 



sin(77r) 
From here we get: 

sin(77r) m 



{mxey + — ^mAe. 



e(A) 



7r(7 - 1) (mA)T 
The expression for p is finally 



A 



m7 



1/(7-1) 



A 



7 



1 



1/(7-1) 



(A - A.)^/^^-^' 



(24) 



(25) 



(26) 



That is, we obtain a power-law behavior with exponent 
/3 = 1/(7 — 1) >1, but now we observe the presence of a 
nonzero threshold 



Ac 



7 ■ 



1 



7717 



(27) 



In this critical threshold reappears in the model. 

However, the epidemic threshold is approached smoothly 
without any sign of the singular behavior associated to 
critical point. 

(iii) 7 > 2: The relevant terms in the Q expansion are 
now 



0(A) 



7 



-mXQ 



7 



7-1 7-2 
and the relevant expression for Q is 



0(A) 



7 



2 1 



7 — 1 A^m 
which yields the behavior for p 
p ^ X — Xc 



X- 



7-1 



(28) 



(29) 



(30) 



with the same threshold Ac as in Eq. ( |27| ) and an expo- 
nent /3 = 1. In other words, we recover the usual criti- 
cal framework in networks with connectivity distribution 



that decays faster than k to the fourth power. Obvi- 
ously, an exponentially bounded network is included in 
this last case, recovering the results obtained with the 
homogeneous approximation of Sec. III. 

In summary, for all SF networks with < 7 < 1, we 
recover the absence of an epidemic threshold and critical 
behavior; i.e. p = only if A = 0, and p has a vanishing 
slope when A — > 0. In the interval 1 < 7 < 2, an epidemic 
threshold reappears (p if A ^ Ac), but it is also ap- 
proached with vanishing slope; i.e. no singular behavior. 
Eventually, for 7 > 2 the usual MF critical behavior is 
recovered and the SF network is indistinguishable from 
an exponential network. 



VI. CONCLUSIONS 

The emerging picture for disease spreading in com- 
plex networks emphasizes the role of topology in epi- 
demic modeling. In particular, the absence of epidemic 
threshold and critical behavior in a wide range of SF net- 
works provides an unexpected result that changes radi- 
cally many standard conclusions on epidemic spreading. 
Our results indicate that infections can proliferate on 
these networks whatever spreading rates they may have. 
These very bad news are, however, balanced by the expo- 
nentially small prevalence for a wide range of spreading 
rates (A <C 1). This point appears to be particularly rel- 
evant in the case of technological networks such as the 
Internet and the world-wide-web that show a SF connec- 
tivity with exponents 7 ~ 2.5 [^,^. For instance, the 
present picture fits perfectly with the observation from 
real data of computer virus spreading, and could solve 
the longstanding problem of the generalized low preva- 
lence of computer viruses without assuming any global 
tuning of the spreading rates jl^,^ . The peculiar prop- 
erties of SF networks also open the path to many other 
questions concerning the effect of immunity and other 
modifications of epidemic models. As well, the critical 
properties of many nonequilibrium systems could be af- 
fected by the topology of SF networks. Given the wide 
context in which SF networks appears, the results ob- 
tained here could have intriguing implications in many 
biology and social systems. 
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